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ABSTRACT 


This thesis presents the theoretical modeling of a terahertz (THz) sensor based on 
a flexible bi-material cantilever that performs thermal deformation under THz 
illumination. This thermal deformation was induced by the large mismatch thermal 
expansion coefficient of two materials. As THz radiation illuminated on the absorber 
plate, the temperature of the structure increased, causing the cantilever to deflect. The 
magnitude of the deflection can be measured electronically or optically, which gives a 
measurement of the input THz signal. The bi-material structure in this thesis is composed 
of a silicon dioxide cantilever beam, with regions coated with aluminum. The effect of 
thickness of two layers on the bending, response time, signal-to-noise ratio will be 
discussed. Maximum bending angle can be found in the results from the example of 
silicon dioxide with thickness layer from 70Onm to 1400nm. The imbalanced stress due 
to fabrication process was taken into account to compare with simulation results to 


provide a better assumption on thickness selection. 
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I. INTRODUCTION 


A. BACKGROUND 


In recent years, research has been rapidly developing in the terahertz (THz) 
frequency range. Due to its non-invasive properties, terahertz radiation has great potential 
in several areas of application such as inspection of fragile artwork, identification of 


concealed weapons, quality control in manufacturing, and medical imaging. 


Detection of THz radiation has been investigated and reported in scientific 
literature [1]. The sensing of THz radiation by the photonic detector techniques normally 
requires a cryogenic operating temperature to reduce thermal noise. The imaging system 
includes a vacuum jacket to produce low pressure and a cooling system, making it 


relatively expensive. 


To overcome this problem, the potential use of un-cooled infrared 
microbolometer camera technology for the THz detection has been investigated [2]. The 
camera’s original focusing optics were replaced with a Tsurupica-based lens and the 
milliwatt-scale quantum cascade laser (QCL) was used as an illuminator. With this 
change, the camera was able to image frequencies down to 0.7 THz. One of the 
drawbacks of microbolometer technology is the generation of heat (self-heating) during 
readout, which limits the integration time to improve the signal-to-noise ratio [2]. In 
addition, microbolometer based focal plane arrays require integrated electronics for 


readout, making them relatively difficult to fabricate. 


An alternative approach that eliminates self-heating as well as integrated readout 
electronics 1s detection based on the deflection of bi-material cantilever structures in 
response to thermal radiation [3, 4]. The amount of deflection can be readout using an 
external light source. These sensors can be easily fabricated using micro- 
electromechanical systems (MEMS) fabrication technologies. The  bi-material 
micromechanical sensor under consideration is based upon a deformed micro-structure 


cantilever under THz illumination. The deformation is caused by different rates of 


expansion due to the difference in thermal expansion coefficients of the two materials. 
Compared to the photon detectors, the bi-material un-cooled thermal detector has lower 
cost and is much more portable. However, the relatively long thermal time constant 
associated with high thermal isolation needed slow results in long response times. It is 
estimated that for real time imaging (30 frames per second), the thermal time constant 
needs to be less than 15 ms. The main focus of this thesis is to optimize the layer 
thickness of thin films used for fabricating the cantilever structure to achieve high 


sensitivity. 


B. BI-MATERIAL CANTILEVER SENSOR 


The bi-material cantilever deforms by absorbing THz radiation by the absorption 
plate (see Figure 1), thereby causing a temperature increase in both materials. Since the 
top material has a higher thermal expansion coefficient than the bottom material, the 
cantilever tends to deform as depicted in Figure | [3]. The magnitude of the deflection is 


then measured electrically or optically and gives a measurement of the input signal. 


Top layer 






Bottom layer-absorption plate 


Figure 1. Deformation of bi-material cantilever (After: [3]) 


The readout system can be either electronic (capacitive or resistive), or optical. 
Electronic readout has been the most commonly used approach, mainly in 
microbolometers. However, as mentioned earlier it has several disadvantages. The main 
problem is its limited integration time due to self-heating. The electronic readout system 


with a large number of pixels requires complex circuitry that will make the 
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microfabrication process very complex and expensive. It will also make it very difficult 
to scale the arrays to high resolution. Compared to the electrical readout system, the 
optical readout is completely external and does not need any electrical leads attached to 
the sensor. This allows the use of low thermal conductivity material, which provides 
higher thermal isolation for achieving greater sensitivity. With this advantage, the optical 


readout becomes more desirable for bi-material based focal plane arrays (FPAs). 


Figure 2 shows the schematic diagram of an optical readout system [4]. The 
incident THz radiation from the target source illuminates the backside of the FPAs 
housed in a vacuum chamber. The absorbed radiation creates an increase in temperature 
resulting in a subsequent deformation of the micro-cantilever pixels in the array. On the 
other hand, a visible light source from an LED illuminates the front side of the FPA. The 


light reflected from the FPA creates a thermal image and 1s recorded by a CCD camera. 


THz source 





THz Lens 











Blade @, 
Wrapped | Beam 
In Plastic | Splitter 
Vacuum | 
Cell 
Aperture 
) Lens1 
LED 
Source 
Imaging 
Screen 
Figure 2. Schematic diagram of an optical readout system (From: [4]) 


As previously discussed, the main physical property of the bi-material cantilever 
is a large difference in thermal expansion coefficients to obtain maximum deflection and 


achieve the operational goal. Clearly, the absorbing material needs to have the proper 


absorption coefficient to ensure the desired wavelengths are captured. The physical 
properties of some of the most commonly used materials available in MEMS technology 


for bi-material thermal detectors are shown in Table 1 [3, 5, 6]. 


Table 1. | Properties of MEMS compatible materials suitable for bi-material sensor 
fabrication 
Thermal expansion Thermal conductivity Heat Emissivity 
coefficient e(Wn'K"!) capacity 6 
a(x10°K™) c(J kg 'K') 
SiNx pie | 19 691 0.8 
S104 0.4 1.4 730 0.8 
Au 14.2 296 129 0.02 
Al 25 231 900 0.02 


As listed in Table 1, SiO. and Al have the lowest and highest thermal expansion 
coefficients, respectively. Comparing physical properties, these materials provide the 
largest difference in thermal expansion coefficients and hence provide maximum 
deflection of the cantilevers. Having a low thermal conductivity, SiO. provides the 
highest thermal isolation, which makes it perfect for use as the structural material for the 


Sensor. 


In the past, there have been difficulties in using SiO. as a MEMS structural layer 
due to the lack of a proper sacrificial layer in the micro-fabrication process [3]. As an 
alternative, SINx was used in combination with Au or Al and became a more popular 


choice of materials for bi-material cantilever structures [3]. 


The S102 based thermal detectors can now be produced using a new fabrication 
process without the use of a sacrificial layer, as reported by Dong et al [7] and as 
discussed in a previous thesis work [3]. Compared to SiNx, overall deflection 
improvement using SiO, was predicted to be 70% [3]. However, the high stress of S102 
creates an extra curvature that affects the reflected LED light signal. This additional 


curvature will remain a challenge but will be discussed in more detail in Chapter IV. 
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I. SENSOR DESIGN AND MODELING 


The structure of the micro-cantilever in this thesis was based on a previous design 
that has been widely used for infrared detectors [8]. The two symmetrical pairs of folded 
legs hold the absorber, which also serves as the reflector for the readout light beam, and 
is located at the center of the detector as shown in Figure 3 [6]. Such configuration 
facilitates the magnification of individual deformations. The total deflection angle of the 
reflector can be calculated as the sum of the individual angles produced on the folded 


legs,ie. 0=6,+6,. 


Deformation 






Supporting frame a sq \, structure 


Top: reflector ¥ YN 





CAI Legs 
i Si02 Bottom: absorber 
Figure 3. The geometry of micro-cantilever with deflection amplification (After: [6]) 


In order to avoid diffraction effects, the size of the individual detectors (pixel) has 
to be larger than the wavelength of incoming radiation that is to be detected. In this thesis 


design, the measurement of one pixel is 200umx200um , which is capable of diffracting 
limited detection incident light from 2 THz (150m) and upward. Another consideration 


of the structure design was the problems associated with partial device release due to the 
imperfections of the etching process. The release involves removal of the substrate from 
the back side using the Bosch etching process [9]. The Bosch process consists of two 
steps: polymer deposition and reactive etching (RIE). As shown below in Figure 4 [9], 


during the RIE step, silicon was removed rapidly by the ion bombardment gas. To 


prevent extensive undercut during RIE (see Figure 4b), repeated polymer deposition steps 


are included to protect the sidewalls as depicted in Figure 4. 


"POLYMER |/_ 





(d) (e) (f) 


Figure 4. The steps of Bosch process where the undercut is minimized using a thin 
polymer layer (From: [9]) 


This process can result in discrepancies between the shape of the mask of the 
opening on one end of the wafer (see Figure 5a) and the actual opening of the mask after 
it reaches the opposite end, as illustrated in Figure 5b. Therefore, the design had to be 


adjusted to account for imperfections in the process. 





(a) (b) 


Figure 5. The over and under etch of different sections during substrate removal to 
isolate pixel 


To overcome this problem, the pixel geometry had to be redesigned while taking 
into consideration the benefits of over etching. In this thesis, instead of having the same 
length of legs in a square area as in Figure 6a, a shorter and outward pair of legs outside 
the pixel was used to preserve extra space for over cutting as shown in Figure 6b. Figure 


7 shows the detailed dimensions of the detector designed by Dr. Dragoslav Grbovic. 





Figure 6. (a) Top view of initial design and (b) top view of new design to reduce over 
etching complete release 





Figure 7. Dimensions of the new pixel design 


A. PERFORMANCE ANALYSIS 


The goal of this research is to improve the performance of the bi-material pixel. 
The performance is described by several parameters. The first parameter is the per-Kelvin 
deflection angle of cantilever where larger deflection angle corresponds to higher 
sensitivity. Therefore, a better performance can be expected due to the signal’s being 
stronger than the surrounding noise. The second parameter is the thermal time constant. 
This describes how fast the detector is. The heated material must cool down to have the 
cantilever return to its original position and be ready for the next operational cycle. This 
process will take some amount of time, which will limit the speed at which the detector 
can be operated. The smaller the relaxation time, the faster detector operation can be. 
Therefore, the goal of this research is to maximize the per-Kelvin deflection angle and 


minimize the thermal time constant (response time). 


To maximize the bending angle due to the given incoming power, the pixel needs 


to have a large increase in temperature, AT, as given by the following equation [3]: 


Tt = (2.1) 


G 


Where P is total absorbed power (watts) and G is thermal conductance (W -K’'), 
which is the ability of material to conduct heat. For the same amount of absorbed power, 
the lower the value of G (value and based on the material and geometry), the higher the 


temperature change and hence larger the deflection angle. 


Thermal conductance also has a great effect on the heat conduction. It is inversely 
proportional to the thermal response time. For the selected material with a heat capacity 
C (J-K'), the higher the thermal conductance G, the smaller (faster) the device’s heat 


dispersion time. The response time 7 1s determined by [3]: 
T=— (2.2) 


As discussed above, the performance of the bi-material device is highly dependent 
on the thermal conductance of the device. The thermal conductance depends on the 
structure geometry design and the materials used. The major challenge in the design of 
the detector is the simultaneous optimization of both the bending angle and the response 
time. In order to maintain simplicity of processing, previous cantilever based infrared 
sensor research [3] provided a starting point for the present work. The detailed 


dimensions of the structure are listed in Table 2. 


Table 2. Structure dimensions used in the analysis 


Pixel size(uum*) Absorber size (um*) Widthoflegs(ym) length of legs (yum ) 
200 x 200 88 x 200 8 see Figure 7 


B. THERMAL CONDUCTANCE AND ITS LIMIT 


The equation for the total thermal conductance for a bi-material thermal detector 


is given by [10]: 


G=G,,+G 


leg rad 


+(G +G 


air ,conv air,cond ) 


(2.3) 


Since the device is typically operated under low pressure in a vacuum cell, the 


convective (G 


air ,conv 


) and conductive (G,,,.,.q ) thermal conductance due to air are 


negligible. The total thermal conductance can therefore be expressed as: 


G=G 


leg 


+G,, (2.4) 


Gieg 18 the thermal conductance through the legs to the heat sink and Graqg is the 
thermal conductance due to radiation emitted from the surface of the pixel. In general, 


Graq 18 relatively small compared to Gjeg The value of Gag is given by [11]: 
G4q = 4AOT (Ey, + Esio, ) Q>) 


where A is the total surface area of the detector, o 1s the Stephan-Boltzmann constant, T 


is the working temperature of the detector, and ¢,, and €,, are the emissivities of Al 


and S10», respectively. 


The general equation for Gjeg is given by [3]: 


G g-W-t 


leg =o ] 





(2.6) 
where g is the thermal conductivity of the material , W is the width of the leg, ¢ 1s the 
total thickness of the leg and / 1s the total length of the legs. 


The equation for the legs composed of Al and SiO» layers is thus given by: 


_ 8a Wty , 8sio, W “tsi, 


Gari sio, - 


(2.7) 


Larssio, Laissio, 
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Since the g,, 1s much bigger than g,,. , the second term can be ignored. The 


equation becomes: 


~ §Al W sty, 
Gar sio, ae (2.8) 
Al+SiO, 


The equation for the legs with only SiQ> layer is given by: 


§ si 2 W obs; 5) 
Gsio, — (2.9) 


lio, 





The G,,, for combining two parallel components in 2.8 and 2.9 can be expressed 
as: 
a ee (2.10) 
Grice Gui. sio, Gio, 


The total G 


eg {Or Symmetric legs on both sides will be: 


t we Tee 
= DW (241 lB sion'Si&, iy (2.11) 


Larssio, § sio,“sio, + Lsio, Sabai 


G 


leg 


Again, since the g,, is much bigger than g,,,., ignoring the 4), so 8sio,tsio, term, 


Equation 2.11 becomes: 


W8 ots, 
eg 2X (2.12) 
Lio, 
Substituting the values of gs. =14Wm"K™), W=8um, Iso = 282um 


from Table 1 and Table 2 into Equation 2.7 will result in G,,, =3.97x10~t,,, (WK '). 


leg 


From Equation 2.5, the Graq can be calculated by using A=1.76x10°(m’), 
assuming the thermal radiation is only emitted from the absorber while the surface area of 
the legs is too small to be negligible, o =5.67x10°(W °K“), T=300K, ¢,, and ¢ 


SiO, 


is 0.02 and 0.08 respectively gives G,, =8.84x10° (WK). 


ra 
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Having a fixed value for the width and the length of the legs, Gjeg can be changed 
by varying the leg thickness. The value of G,ag depends only on the surface area of the 
pixel, which is nearly independent of the thickness of the legs. Thus, for a fixed pixel 
area, the G,aq remains constant while the Gjeg decreases with its thickness. Therefore, the 
minimum value of thickness to reduce the overall thermal conductance can be obtained 


When Gijeg is equal to Gyaq as illustrated in Figure 8. 


However, due to the issues of stress differential resulting from the micofabrication 


process, 700 4m is the lowest value of thickness of SiO, that has been used in these 





calculations. 
@ 
o 
c 
is] 
—_ 
oO 
5 
To 
= 
© 
oO 
rs 
= 
= 
rT) 
= 
- 
0 tnin= 2.22 um 
Thickness 
Figure 8. Thermal conductance as a function of leg thickness 
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HI. FINITE ELEMENT SIMULATION 


The response of the pixel for absorbed THz power is simulated using COMSOL 
finite element simulator. In the model, the top reflector material was temporarily removed 
in order to eliminate the stress effects in the central region of the pixel. In addition, the 
exact material to be used for both absorption of THz radiation and reflection of the 
readout is a subject of a separate research project. The value of incoming power in this 
simulation was estimated based on a previous experiment where | mW incident power 


was illuminated on 120x160 pixel infrared FPA with 504mx50um per pixel area. This 


corresponds to an intensity of incident power of about 20.8Wm~ and in COMSOL 


simulation, 100% absorption is assumed. 


It is important to determine the resonance frequencies of the entire structure since 
it can affect the speed of operation of the pixels under time varying THz power. If the 
lowest resonance frequency is close to the frequency of the time varying THz signal, a 
non-linear response can result from the sensor. On the other hand, higher resonance 
frequencies are associated with stiffer spring constant resulting in smaller deflection of 
the cantilever. Figure 9 shows the first three eigen modes simulated using COMSOL of 
the pixel structure shown in Figure 7. The first mode is the most important in the design 


of the pixel. 





(a) (b) (C) 


Figure 9. The lowest three eigen modes (a) 4921 Hz, (b) 5707 Hz and (c) 15104 Hz 
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The frequency response of a resonance is usually not be exactly symmetric about 
the resonant frequency. The amplitude of the oscillations x for a damped oscillator can be 
expressed as: 

Fim 


aa (3.1) 
(a, -@ YP +(yoy 


Where F is the driving force, m is the mass, @ is the driving frequency, @, 1s the 
resonant frequency and y is the ratio of damping coefficient b to the mass m. A typical 


resonant curve is shown in Figure 10 (a). When the resonant frequency is much larger 


than the driving frequency (1.e., @< q, ), the equation can be simplified as: 


ver im_Fim_F 
oO; kim ok 





(3.2) 


which shows that the response amplitude is proportional to the driving force, F =kx, 
where k is the spring constant. Under this condition, the amplitude of oscillations is 
independent of the driving frequency and proportional to the driving force providing a 


linear response. 


In the design, the first mode for different leg thicknesses was simulated and 
graphed in Figure 10(b). It was found that the lowest resonant frequency was 2000 Hz. 
This result is far from the speed of operation of the pixels (30 Hz) needed for real-time 
imaging. It is known that a frequency of 30 Hz is the minimum speed for a human eye to 
decipher a motion picture. The thickness in the simulation will present a steady condition 


for the amplitude response. 
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Figure 10. (a) Amplitude as a function of frequency for a damped oscillator and (b) 
resonant frequency of a pixel as a function of leg thickness 
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A. MAXIMUM BENDING 


As shown in Figure 11, the total displacement of the cantilever in the z-direction 
can be calculated by summing the displacement at points | and 2. The sine of the 
deflection angle (@) is trigonometrically equal to the total displacement in the vertical 
direction divided by the length of the pixel. For small@, sin@ ~ @ and the deflection 


angle are equal to the vertical displacement. 


Point 1 
ia 





Point 2 


Figure 11. | Determination of angle of deflection 


By measuring the temperature increase in the absorber (the temperature difference 
. d . 
between heat sink and absorber), “ can be calculated for different leg thicknesses. 


Aluminum has a higher thermal expansion coefficient than silicon dioxide and, thus, 
increasing the thickness of the Al will increase the bending angle as shown in Figure 12 
(a). However, it was found that for Al thickness above 400 nm (see Figure 12 (a)), the 
bending angle starts to decrease. This is primarily due to reduction in temperature rise of 
the pixel as the aluminum thickness is increased as a result of higher thermal conductance 


of the legs. The smaller temperature increase incurs a smaller deflection. The value of 


< then can be simply obtained by dividing one by the other, as shown in Figure 12 (b). 
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The - curve is almost identical to the bending angle curve. As the author fixed the 


thickness of Al and varied S102, the = curve had the same behavior. 
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Figure 12. = (a) Deflection angle and change in pixel temperature as a function of Al 
thickness with 900 nm SiO», (b) dependence of “ as a function of Al thickness with 


900 nm Si0O> 
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The 3D graph in Figure 13 shows the tendency of the ae curve by varying both 


aluminum and silicon dioxide thicknesses. The data in Figure 13 can be used for 


determining the optimum parameters for achieving the highest rate of deflection with 
temperature. 





0.006 
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¥ 0.004 
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iS 0.003 
S 
0.002 
Figure 13. 


3D graph of the < by varying the thicknesses of both SiO) and Al 
B. RESPONSE TIME 


In order to evaluate the time constantz , Equations 2.2 and 2.6 can be used along 
with the heat capacity of the absorber (200 x 88 um) and thermal conductance. By 
assuming the mass of the Al layer is too small to be negligible, the heat capacity of the 
structure can be calculated by the specific heat capacity of SiOz multiplied by the mass of 
SiO, obtained from COMSOL, C (JK ') = 730 (JKg 'K ')xmass (Kg), which means the 
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heat capacity is directly proportional to the total thickness of the S10, layer since the 
pixel area is a constant. As mentioned in Chapter II, the thermal conductance is 
proportional to the thickness of SiO and Al. As a result, by fixing the thickness of the 
S102 layer and increasing the thickness of the Al layer, the time constantz decreased as 
shown in Figure 14. However, fixing the thickness of the Al layer and increasing the 


thickness of the S102 layer, the time constantz increased. This indicates that the heat 


capacity dominates the thermal conductance in optimizing the response time. 





Figure 14. | 3D graph of the response time by varying the Al and SiO» thicknesses 


C. QUALITY FACTOR 


The quality factor, Q is necessary in the calculation of the Noise Equivalent 


Temperature Difference (NETD). The details of the NETD and its computation will be 


discussed in the next section. Q-factor of a resonator is generally defined as the ratio of 
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the maximum energy stored to the energy lost in a single oscillation. Therefore, devices 
with higher quality factors will have larger stored energy when driven at the resonant 


frequency. 


Q factor also represents the sharpness of the resonance curve shown in Figure 16. 


A high quality factor is characterized with a high response peak and a narrow bandwidth, 


. 1 
which is the frequency difference between the points of amplitude at —— (also known as 


J2 


3dB) of the peak, as shown in Figure 15 [12]. 
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Figure 15. | Bandwidth of resonant frequency (From [12]) 


Quality factor can be expressed as: 


_ energy stored percycle — ——_—~ resonant frequency 


energy dissipated per cycle handwiktiut i of peak 


J2 


(3.3) 


The damping constant (Cgamping) due to the motion of a cantilever in air can be 


estimated using [13]: 
Camping = Cantilever width -1./bpap (3.4) 
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where wis the vibration frequency, p is the fluid density, and yw is the fluid viscosity, 
b is a constant number dependant on the structure design. In air, the density 


is 0 = 1.025 kg-m°~ and viscosity is u=1.871x10°Pa:s. 


Equation 3.4 was based on a cantilever with rectangular cross-section moving in 
the air where the parameter bis approximately equal to 2. Considering the multi-leg 
design of the structure, a higher damping and hence a higher value of b is expected. The 
increase in b is due to the increase in the perimeter of the structure resulting in increased 
air damping as it vibrates. The value of b is usually estimated by measuring the 
bandwidth of the resonance peak experimentally. Figure 16 shows the experimental 
setup used for measuring amplitude as a function of driving frequency using a laser 
vibrometer. A speaker was used as a sound source, providing an acoustic plane wave in 
the 100 to 7000 Hz frequency range. In the measurement, a single pixel of a FPAs 
fabricated is employed with laser light from the vibrometer, focused using a microscope 


objective lens assembly. 





Figure 16. | Lab setup for measuring frequency response of a pixel. (A) Speaker for 
sound source with directional cone. (B) Laser vibrometer with a microscope objective 
lens. (C) FPAs 
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As seen in Figure 16, the measured frequency of a pixel is compared with the 


simulation with b = 3 and without residual stress (green line), showing a higher resonant 
frequency and a symmetric peak. However, a residual stress of 10’Pa on Al and 


3x10°Pa on SiO» was found to exist by the characterization of thin films performed 
during the microfabrication process. Inclusion of the stresses in the simulation resulted in 
decrease of the resonant frequency; the curve became asymmetric (blue line in Figure 17), 
which is closer to the experimental data but with a relatively wide bandwidth. The red 
line in Figure 17 is the simulated data with b = 2 and stresses that show a close agreement 


with the measurements. 


—— Experimental data 
b=3, Al= 200 nm without stress 
»=— b=3, Al = 200 nm with stress 
*— b=2, Al= 200 nm with stress 


Normalized Amplitude 





2000 3000 4000 5000 6000 
Frequency (HZ) 


Figure 17. | Comparison of simulated frequency response with experimental data 


By using Equation 3.4 with b = 2, the frequency response is simulated for a set of 
pixels with different layer thickness parameters; the results are shown in Figure 18. The 


value of the Q factor can be calculated by Equation 3.3 and ranges from 20 to 45. These 


Z2 


values are then used to calculate the Noise Equivalent Temperature Difference (NETD) 
performance of the detector in the FPAs for the next section. Note that the Q factor in 


operating condition (vacuum) will be higher. 
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Figure 18. | Frequency responses of pixels with Al thickness from 160 nm to 200 nm 
and S10, thickness from 800 nm to 1400 nm 
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D. NOISE EQUIVALENT TEMPERATURE DIFFERENCE 


Noise Equivalent Temperature Difference (NETD) of a sensor is defined as the 
amount of temperature of a target above the background needed to produce a signal-to- 
noise ratio equal to | at the sensor. Both the target and background are assumed to be 


blackbody sources. An equation of NETD can be expressed as [3]. 


NETD, =< (3.5) 


AZ 


where Zn is the amplitude of noise and AZ is the amplitude of signal. For the micro- 
cantilever, the signal will be the deflection of the device as a result of power from the 
target, which can be expressed as [3]: 


a5 T (3.6) 


GFE OO pant 


where A is the surface area of a pixel, F is the f-number of the lens used for focusing the 


radiation beam onto the detector, and 7, is the transmission of the lens. In the case of THz 


detection, the lens is made of polyethylene. The average transmission of the lens in the 


THz region (1-5 THz) is about 0.77, as shown in Figure 19 [14]. 
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Figure 19. = Transmission of 2 mm-thick high density polyethylene in THz region (From 
[14]) 


. = dP 
R(w) 1s the responsivity of the sensor, T is the room temperature and CP 4-4, 18 


the slope of the black body radiation curve from /, to 4, . For the desired frequency range, 


1-5 THz, & ,,-4, can be obtained by approximating the power emitted by a black body 


source in THz spectral range as [2]: 


P(T) = 


V 
20kT 
[av 2S ve (3.7) 
C 


0 


where v, is the upper limit of frequency range and the lower limit is assumed to be 0. By 


taking the derivative with respect to T on both sides, the equation becomes: 





= aera et (3.8) 
C 


Using cutoff frequency v, = 5THz and calculating & from Equation 3.8 gives 


a value of 0.06. 


The primary source of micro-cantilever noise comes from the temperature 


fluctuations. The root mean square (average) of this noise can be expressed as [3]: 


R(o@)T 4k, BG 
aS ie a aE nde na (3.9) 
7] 


Where k, is the Boltzmann constant of 1.38x10°m*kgs °K” , B is the frequency 


bandwidth of the sensor—which is limited by the thermal time constant (~ 30Hz )—G is 


the thermal conductance of the sensor and 77 1s the absorbance (~0.3) for a typical THz 


absorbing material used in sensor fabrication. 


Additional noise is introduced from the lattice vibrations due to the thermal 
energy (kg7), which excite some of the resonant modes. This additional vibration is called 


thermo-mechanical noise and is expressed as [3, 15]: 


(3.10) 





Notice that the operating frequency @ in this case is much smaller than the 


resonant frequency @ ; therefore Equation (3.10) can be simplified as: 


<Ziy > = (3.11) 
Mp 


where Q is the quality factor and k is the spring constant. 
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By the definition of NETD (Equation 3.5), the NETD,, (due to temperature 
fluctuations) and the NETD,,, (due to thermal mechanical noise) then can be obtained 


by combining Equation (3.9) with (3.10) and Equation (3.9) with (3.11): 





AF°T./4K,,BG 
NETD,,, = — (3.12) 
TIA a~m 
2 
NETD,,, -—* = ae (3.13) 
TyAR(@\( Tan, VON 


The total value of NETD thus can be evaluated by taking a square root of the sum 
of the squares of individual NETDs. 


The graphs of NETD with - and NETD with response time are shown in 


Figures 20 (a) and (b), respectively. As displayed in the graph, when the thickness of Al 
is increased, the value of NETD increases slowly at about 0.02 mK per 50 nm. When the 
thickness of SiO, increased, the value of NETD increased rapidly at about 2 mK per 50 


nm. 
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Figure 20. (a) NETD with = (b) NETD with response time 
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IV. ANALYSIS OF RESULTS 


A. SELECTION OF THICKNESS 


It 1s important to select the most favorable thickness combination of the bi- 
material detector in order to achieve the optimal performance. The challenging part 
comes from the need to simultaneously optimize all the desirable parameters such as 
sensitivity, response time and noise. As shown in the equations above, these complex 
physical properties have adverse effects on each other. Therefore, optimizing one often 


worsens another. 


One way to simplify this problem is to create a simulation of numerous values of 
different thickness settings and then graph the results. The desired properties of 
unsimulated data can be predicted from the resultant curve. From previous chapters, the 


3D graph in Figure 14 shows that the response time will decrease by either reducing the 
thickness of SiO2 or increasing the thickness of Al. The maximum value of = can be 


found in each curve from the 3D graph on Figure 13, which was simulated by 100 nm 
intervals of the SiO, layer thickness. The peak point with the highest deflection 
amplitude and the smallest total thickness with the smallest response time from Figure 21 
will provide the best solution. In conjunction with these selection criteria, the least total 
thickness that can be considered must be 700 nm due to restrictions on current fabrication 
techniques. As discussed previously, the value of NETD of all the combinations with 
S102 thicknesses less than 1400nm will never be higher than 62 mk, which is relatively 


low for uncooled thermal sensors in the THz spectrum. 
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Figure 21. | Response time and peak deflection rate as a function of SiOz thickness 


Another method for determining the combination of thicknesses is to fix one or 
two variables and vary the remaining parameters. For example, in Figure 22, if the 
desired NETD must be lower than that of the thickness combination of 4 (SiO 800 nm, 
Al 200 nm), which is 46 mK, and the response time has to be lower than thickness 
combination | (SiO. 800nm, Al 80nm), which has a value of 0.75 second, this will limit 
the choice of bending angles of combinations 2 and 3. The combination of 3 will thus be 


selected because it provides the maximum angle. 
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Figure 22. | Thickness sets with the optimized desirable parameters 


B. THE EFFECT OF STRESS FORCE 


Another important physical effect on the performance of bi-material cantilevers is 
the initial deflection due to the residual stress in the two materials. This stress builds up 
during the chemical deposition stage of the micro-fabrication process [16, 17]. The 
residual stress is a result of a difference in the shrinking rate between two materials in the 
x-direction and y-direction when cooled down. These effects form a curvature as shown 


in Figure 23. 
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Figure 23. | Top and side views of (a) structure without residual stress and (b) structure 
with residual stress 


The initial curvature highly affects the optical readout caused by thermal 
deflection due to THz absorption. These factors are important for consideration in 


performance analysis in order to properly choose the realistic thicknesses for simulation. 
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V. CONCLUSION 


This thesis has successfully demonstrated the use of COMSOL multi physics 
finite element simulation with heat transfer and solid stress models for the design of a bi- 
material cantilever sensor in the THz region. The results for changing thicknesses and the 
physical effect on the sensitivity, response time and noise were reported. From these 
results, the best combination will be compared experimentally to the proposed device. 
The characterization methods developed in this study laid a foundational process for a 
more advanced structure design and more complex simulations. In the future, coupling 
these results with the COMSOL RF model could enable a simulated heat generation 


source from THz radiation to more closely emulate the actual device. 
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